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Abstract. The large-scale gradient of the interstellar 
medium (ISM) density distribution essentially affects the 
evolution of Supernova remnants (SNRs). In a non- 
uniform ISM, the shape of SNR becomes essentially non- 
spherical, and distributions of gas parameters inside the 
remnant become strongly anisotropic. The well-known 
self-similar Sedov solutions may not be applied to mod- 
elling such non-spherical objects. Therefore we propose a 
new approximate analytical method for full hydrodynam- 
ical description of 3D point-like explosions in non-uniform 
media with arbitrary density distribution. 

On the basis of this method, we investigate the general 
properties of evolution of 2D non-spherical adiabatic SNRs 
in ISM with large-scale density gradient. 

It is shown that the real shape of adiabatic SNR be- 
comes more non-spherical with age, but the visible shape 
remains close to spherical even for strong real asymme- 
try (ratio of maximal to minimal shock radii) and surface 
brightness contrast. It is shown also that the values of the 
parameters of X-ray radiation from the entire SNR (lumi- 
nosity, spectral index) are close to those in the Sedov case 
with the same initial data. However, the surface distri- 
bution of the X-ray emission parameters is very sensitive 
to the initial density distribution around the SN progen- 
itor. Therefore, the X-ray maps give important informa- 
tion about physical conditions inside and outside of the 
non-spherical SNR. 

Key words: hydrodynamics - method: analytical - ISM: 
supernova remnants - X-rays: interstellar 



1. Introduction 

The investigation of Supernova remnants (SNRs) gives im- 
portant information on the physics of Supernova (SN) ex- 
plosions, on the properties of the surrounding interstellar 
medium (ISM) and on shock wave physics. The majority 
of galactic S NRs are in the adiabatic stage of evolution 
(Lozinskaya 1992 ). If the density of the ISM is uniform, 
their hydrodynamics are well described by the self-similar 
Sedov solution (Sedov |1959|, Shklovskiy |l962[). The typi- 



cal values of plasma temperatures in SNRs are T ~ 10 to 
10 s K and, therefore, SNRs radiate mainly in the X-rays. 
The spectral characteristics of the equilibrium X-ray emis- 
sion for a plasma typical SNR abundances was calculated 
by Shapiro & Moore ( |1976| ), Raym ond fc Smith ( |1977D , 
Shull ( |l98ll ) and Gaets & Salpeter ( |l983| ). 

But the real situation is more complicated. In many 
SNRs, the plasma is in noncquilibrium ionization (NEI). 
Often there is no thermal equilibrium between the elec- 
trons and the ions. The physical conditions in the inner 
parts of SNR may be modified by electron thermal conduc- 



tivity (Itoh 1977, Cox fc Anderson 1982, Hamilton et al. 



1983| , Jerius & Teske |1988j , Borkowski et al. |1994| , Bocchino 
ct al. |1997p . There are also other effects which affect the 
plasma emission in SNR, but their influence is small: res- 
onant scattering, diffusion etc. (Raymond & Brickhouse 



1995 and references there). The Sedov solution may also 



be modified by the presence of sm all-scale cloudlets in the 
ISM (Byc hkov & Pikelner |1975|, McKee & Cowie |l975, 



Sgro |1975|, White & Long |l99l|) . 

Practically the all above-mentioned investigations as- 
sumed spherical symmetry of SNRs, as the result of 
their evolution in uniform ISM. However, the observa- 
tions show predominantly nonsphcrical shapes (Seward 
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1990| , W hiteoa k & Green |1996[ ) . For example Kcst even & 
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a separate class of barrel-like SNRs. We have carried out 
the analysis of visual anisotropy for SNRs from the cat- 



alogues of Seward ( [L990[ ) and Whiteoak & Green (|99§) 
(Fig. [I]). This figure shows that typical values of the visual 
anisotropy is usually d max /d mul — 1 to 2. 

Many SNRs with nearly spherical visual shapes have 
an anisotropic distribution of surface brightness (e.g. Ke- 
pler SNR, Cygnus Loop, RCW86 etc.). 

Therefore, truly spherical SNRs are considerably rarer 
than believed. 

Non-spherical SNRs (NSNRs) may be created by a n 
anisotropic SN explosion (e.g. Bisnovatyi-Kogan 1970 ). 
Non-spherical shapes of SNRs in the free expansion stage 
(Fig. |l|) should be mainly produced in this way. Another 
important reason for the non-sphericity of adiabatic SNRs 
can be a non-uniform density distribution of the ISM or 
the large-scale magnetic field. In such cases, self-similar so- 
lutions cannot be used, while direct numerical calculations 
of the problem are difficult because of the complications 
of 3D hydrodynamical modelling SNR evolution multiplies 
by the complications of nonequilibrium describing of gas 
element evolution. Therefore, at present only a few sim- 
plified models have been built (Tenorio-Tagle et al. 1985j , 
Bodenheimer et al. |1984|, Claas & Smith |l989j, Bocchino 
et al 



1997) 



Real possibility to perform an investigation of the evo- 
lution of non-spherical SNRs bases on approximate meth- 
ods for hydrodynamics. The thin-layer (Kompaneets 1960) 
approximation for the calculation of SNR shape is widely 



used (Lozinskaja 1992, Bisnovatyi-Kogan & Silich 1995 



and references there) . But this approximation has low ac- 
curacy for the adiabatic stage in a nonuniform medium 
and does not allow to calculate the behaviour of the gas 



inside th e SN R (Hnatyk |1987| , Hnatyk & Petruk |1996[ ). 
Hnatyk (1987) has shown that for our probleme it is more 
promising to develop approximate methods under a sec- 



tor approximation (Laumbach & Probstein 1969), which 
allows to calculate both the SNR shape and the ga s char- 
acteristics. In the work of Hnatyk & Petruk ( 1996| ) a new 
approximate analytical method for the complete hydro- 
dynamical description of a point explosion in a medium 
with an arbitrary regular but smooth density distribu- 
tion was proposed . It combines the advantages of the two 
above-mentioned methods and, therefore allows to calcu- 
late the hydrodynamical aspects of non-sherical SNR evo- 
lution with high enough accuracy in a short computing 
time. 

We use this method here, but we limit ourselves to the 
case of equilibrium emissivity. Non-equilibrium effects will 
be considered in a further paper. 



2. Hydrodynamical modelling 



As proposed by Hnatyk & Petruk ( 1996 ) the hydrody- 
namical description includes two steps: the calculation of 
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Fig. 1. Distribution of ratios of maximal to minimal visual 
diame ters o f SNRs from 1 - radio observations (Whiteoak & 
Green 1996) and from 2 - an X-ray catalogue (Seward 199C). 



The young SNRs, which are expecting to be in the free e xpan- 
sion stage, ar e sep arated - 3 (Lozi nskaya 



1993, Seward 
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Greiner et al. |l994| , Aslanjan |l996j Stankevich 1l996| ) 
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shock front dynamics (shape of SNR) and the calculation 
of the state of the plasma inside the SNR. 

2.1. Calculation of shock front motion 



Klimishin & Hnatyk ( |l98l[ ), Hnatyk ( J1987| ) have shown 
that the motion of a strong one-dimensional adiabatic 
shock wave in a medium with an arbitrary distribution 
of density is described with high accuracy by the approx- 
imate formula: 



dR 
~dt ~ 

where 



D{R) = const • (p°(R) ■ R N+1 ) " 



, = f 1/2, m{R) <N + 1 
" < 1/5, m(R) > N + 1 



(1) 



(2) 



R is the distance from the explosion, p°{R) is the initial 
density distribution of the surrounding medium; m(R) = 
-d\n(p°(R))/d\nR\ N = 0,1,2 for a plane, cylindrical 
and spherical shock, respectively. 

Formulae ([j])-(||) generalize two basic features of shock 
motion in non-uniform media: 

- a deceleration in a medium with increasing (m{R) < 
0), constant (m(R) = 0) or slowly decreasing density (0 < 
m(R) < N + 1); then the parameter k is close to k = 1/2; 

- an acceleration if the density is decreasing fastly 
enough as m(R) > N + 1; then the parameter k is close 
to k = 1/5. 

Since in all realistic cases m(0) = 0, i.e. the initial stage 
of the motion of the shock from a point explosion is always 
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described by the self-similar Sedov solution for a uniform 
medium, for the deceleration stage of shock motion from 
cq. (0)-(|) we have: 



D(R) 



2 



E n 



3 + A \a A (N,7) ■ (f>(R) 



1/2 



■R 



-{N+l)/2 



(3) 



where a A {N, 7) is the self-similar constant for a uniform 
medium (a A (0, 5/3) = 0.6029, a A (2, 5/3) = 0.4936), E Q 
is the energy of explosion. 

If the density distribution is such that a region of ac- 
celeration with m(R) > N + 1 exists and begins at some 
distance R\ where m(i?i) = N + 1, then an approximate 
formula for shock velocity is: 



Da{R,Ri) ~ Du(Ri) 



p°(R) ■ R N + 1 



N+l \ 1/5 



(4) 



Here, Dd(-Ri) is determined from (||). 

If the density distribution is such that there is a transi- 
tion from a region of acceleration with m(R) > N + 1 into 
a region of deceleration where m(R) < N+l, then another 
feature appears. At the beginning, the shock deceleration 
will be analogous to the motion of an external (forward) 
shock as described by Nadyezhyn (1985) and Chevalier 
( |1982| ). The external shock decelerates with k « 1/5 and, 
when its velocity equals velocity D^{R) from (||) the fur- 
ther deceleration is described by k w 1/2. 

Therefore, a general formula for the shock v elocit y 
which takes into account the three cases is (Hnatyk |l987 ): 



D(R) 



D D (R), 
min[D A (R- R 2n -i); D D (R)}, 



R 



2n 



< R < Ri 

_1 < R < i?2n^ 



(5) 



where n = 1, 2, 3...; R n are zeros of the function m(R) — 
(N+l) (the points of changing of regimes of shock motion) 
in increasing order of R in the interval (0, 00). 

It must be emphasized that the majority of models 
interesting for astrophysics are described more simply by 
formulae ©-§). 

The equation of the trajectory of the shock motion is 



dR 
D(R)- 



(6) 



When the density distribution departs from spherical 
symmetry, the calculation should be carried out by sectors. 
The 3D region is divided into the necessary quantity of 
sectors and Eq. (||) is integrated for each of them. Such 
an approach allows also to consider an anisotropic energy 
release E(Q) replacing E a by 2ttNE(Q) for N = 1, 2. 




Fig. 2. Profiles of a shock front in a non-uniform exponential 
medium Eq. ([ll]) as a function of dimentionless time r (see 
text). The coordinates are non-dimensional as denned in the 
text. Line 1 is the result of the numerical 2D calculation; line 
2 is the approximate law (^); line 3 is the result of numerical 
calculation in the sector approximation; line 4 is the result of 
the thin layer approximation. 



2.2. Calculation of the plasma characteristics inside SNR 

The second step concerns the determination of gas pa- 
rameters inside the SNR. It uses on the shock motion law 
by sector discribed in the previous step. We will use the 
Lagrangian coordinates (a, t) (a is the initial coordinate 
of the gas element at t = 0) and take as the unknown 
function the position (Eulerian coordinates) of the gas el- 
ement r(a, t) (0 < r(a, t) < R) at time t. The function 
r(a,t) is expanded into series about the shock front and 
about the center of explosion. Then these two decomposi- 
tions are combined. The coefficients of decomposition near 
the shock front is determined by the motion law in each 
sector and near the center of e xplo sion assuming a zero 
pressure gradient (see Appendix |A.3| ). It is important that 
the value of pressure in the central region be taken equal in 
all sectors contrary to previous propositions (Laumbach & 
Probstein |l969i Hnatyk [l987l Hnatyk |l98^ ) . The value of 




Fig. 3. Distributions of density (a, b) and pressure (c) inside an SNR expanding in a fiat exponential medium Eq. ( |ll| ) for 
different directions from the place of explosion at the dimensionless time r = 4.268. The solid lines corres pond to the present 
method, the dashed lines are the result of the direct 2D numerical calculations of Kestenboim et al. (1974). Here 7 = 1.4. 



pressure around the point of explosion is determined from 
the condition that the ratio 77 of pressure in the center of 
explosion P(0, t) to the average value of energy density 
inside SNR E /Vt t(t) does not change with time and is 
equal to the self-similar value for a uniform medium: 



V = P(0,t)V to t(t)/E o = rj A . 



(7) 



In this way, we approximatively take into account a redis- 
tribution of energy between the sectors in the cases of an 
anisotropic explosion and/or a non-uniform medium. 
The calculation of r(a,t) is presented in Ap- 



pendixes A.l- A. 3 



The other parameters (pressure P (a, t), density p(a, t), 
velocity u(a, t)) are exactly derived from r(a, t) in case of 
adiabatic motion of gas behind the shock front (Hnatyk 



1987) 



- the density is obtained from the continuity equation 



(8) 



- the pressure is derived from the equation of adiabaticity 
P = Kpi or 



P(a,t) 



p°(a) 



( D{a)V ( p{a,t)^ 
\D(R)J \ P (R,t), 



P(R,t) \p°(R), 
- velocity directly from approximation r — r(a, t) 
dr(a, t) 



it(a, i) 



dt 



(9) 



(10) 



2.3. Testing the method 

The proposed method exactly describes the shock trajec- 
tory in self-similar (Sedov) cases with m{R) — const < 
N + 1, including in particular the uniform (m = 0) 
medium. The space distributions of the gas parameters: 
pressure P(r,t), density p(r,i), temperature T(r,t) and 
velocity u(r, t) inside the SNR are accur ate within 3% in 
the m = case (Hnatyk & Petruk 1996 ). 

The result of the calculation for a point explosion in 
medium with an exponential density distribution 



p°(z)=p°(0)-exp[-- 



z = r cos f 



(11) 



where r is the distance from the center of explosion, 9 is 
the angle between the considered direction and the di- 
rection opposite to the density gradient, and H is the 
scale height, are presented in fig. |^-[| The simulations 
are carried out for dimensionless parameters r = t/t m , 
r* = r/R m , p* = p/p m , P* = P/P m , with t m = 
a A (2, 7 ) 1/2 K~ 1/2 po /2 (0)i?m /2 , R m = H, p m - p°{0), and 
P m = Q!a(2, -f)~ 1 E a R~ i 3 . They allow to obtain the solu- 
tion for any value of the initial parameters E Q , p°(0), and 
H. 

Fig. U presents the results of a direct 2D numerical 
calculation together with result s of d ifferent approximate 
approaches (Kestenboim et al. 1974 ). As we see, for the 
largest time r = 17.9, the proposed method is accurate 
within 7%, i.e., of the same order of accuracy as the nu- 
merical method. Even a numerical calculation by sector 
(ID numerical calculation in each sector) yields a consid- 
erably lower accuracy, only about 20%. 



i l i ~] — i—i— i — r~t—i — i — r~i~^ — I i i i I i i i I — i — i — I — i i i i i i — i i i I — r~m — | i i i | i i i | i i i | i i i | — i — i—r 

-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 

X* X* 

Fig. 4. a) Profiles of shock front expanding in a exponential medium Eq. (solid lines) and in a uniform medium (dashed 
lines) for dimensionless times r = 1, 3, 7. The dotted line is the profile fitting by a sphere for r = 7. b) The 3D shape of the 
same SNR for r = 7. The angle between the symmetry axis of the SNR and the plane of sky is 45°. 



Another important quantity useful to check the accu- 
racy the shock acceleration in exponential atmosphere is 
the breakout time Tb r , when shock velocity D(8 — 0°) — > 
oo. The combination of the 2D numerical solutions for 
t < 17.9 with the self-similar solution for r > 17.9 for 
7 = 1.4 yields Tbr = 29.6 (Kestenboim et al. 1974) while 
approximation (|^) gives Tbr = 30.8. 

The accuracy of the density and pressure calculations 
in the proposed method is illustrated by Fig. ||. 

Additional results of testing of the method for more 
complex density distributions, such as a density disconti- 



nuity are presented by Hnatyk (1987). All these results re- 
veal that the proposed method has high enough accuracy 
in all cases in the sector aproximation. It fails however in 
the case of shock interaction with a small dense cloud. 

3. SNR shapes in non-uniform media 

We now consider the role of the surrounding ISM on the 
evolution and the X-ray emission of adiabatic SNR on the 
examples of media with flat exponential and spherically- 
symmetrical power-law density distributions. 

3.1. Shapes of SNRs in a flat exponential medium 

The exponential law distribution is frequently encountered 
in nature, especially in galactic disks. We consider the evo- 
lution of the shape of the shock front from a point explo- 
sion in a medium with exponential density distribution 
Eq. (0). The morphological evolution of such SNR are 



presented in Fig. [|. We can see from this figure the re- 
markable insensitivity of the visible form of the SNR to 
the ISM density gradient. 

The apparent center of the non-spherical SNR does 
not coincide with real progenitor position. This may be 
important for localizing a possible compact stellar rem- 
nant (pulsar or black hole). 

The evolution of some shock characteristics is pre- 
sented in Fig. ||. The main result is that the average visible 
morphological characteristics of the non-spherical SNRs 
are usually close to those for Sedov SNRs with the same 
initial parameters. 

3.2. Shapes of SNRs in a power-law medium 

Another widely-used density distribution is the power-law 
one, created by stellar winds, previous SN explosions etc.: 



(12) 



We consider the evolution of the shock from a point explo- 
sion in a medium with a spherically-symmetrical power- 
law density distribution when the explosion point is dis- 
placed by a distance r from the center of symmetry f = 
(wind source etc.). Therefore the density distribution as a 
function of the distance r from the explosion point r Q is: 



p°(r,0)=p°(O)l 



2rr„ cos 9 



(13) 
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Fig. 5. Evolution of the shape and characteristics of a shock wave for a point explosion in a flat exponential medium Eq. (|ll|). 
a: 1 - radius of shock front R%/ 2 in direction 6 = 90° (this equals R s for the case a uniform medium); 2 - Rq; 3 - R^; 4 - 
(Ro + R%)/2; 5 - average visual radius of NSNR in exponential medium, b: 1 - ratio of maximal to minimal diameters of visual 
shape of shock wave front; 2 - ratio Ro/R-n', 3 - \og(p 7T / po); 4 - log(To/T x ). 




Fig. 6. a and b The same as in Fig. ^ for a power law density distribution Eq. ( |l3| ) for r ( * = 1.75 and r = 7. The angle between 
the symmetry axis of NSNR and the plane of the sky in the right figure equals 15°. 
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Here p°(0) is the initial density in the point of explosion. 
Hereafter, we take w = —2 and r = 1.75i? m . 

The shape evolution of the SNR in this density distri- 
bution is shown in Fig. ^. It is worthy to note that visible 
shapes of such SNRs may be elongated transverse to the 
density gradient. 

As one can see from Fig. |g, b like to the previous case 
(Fig. [|, b), the projection of the SNR on the plane of the 
sky can cause a spherization of the visible SNR shape. So, 
even a visible spherical shape of SNR does not guarantee 
the uniformity of ISM and isotropy of explosion. 



3.3. Discussion 

From the above results, it follows that: 

1. The non-uniformity of the surrounding medium 
causes asphericity of SNRs. The visible shapes may be 
elongated not only along (Fig. ||) but also transverse to 
(Fig. ^) the density gradient, depending on the type of 
density distribution. 

2. SNR may have an apparent shape close to spherical 
even in cases of essential anisotropy of the real form and 
essential gradient of density distribution along the surface. 

3. The observed anisotropy of the shape is smaller than 
the real anisotropy as result of projection. The visible 
shape remains, generally, close to spherical. 

4. The non-uniformity of the surrounding medium re- 
sults in differences of shock characteristics along the SNR 
surface. If the initial density p varies along the shock sur- 
face, the maximal contrasts are expected in the X-ray sur- 
face brightness (oc p 2 ), they are smaller for the temper- 
ature distribution (oc p) and minimal in shock and post- 
shock gas velocities (oc p -1 ^ 5 ) (Figs. |] ||). 

Therefore to determine the real conditions inside and 
around SNR, it is necessary to use additional information 
about SNR. We consider further the X-ray observations 
as an effective tool for SNR diagnostics. 



4. Integrated characteristics of the X-ray radiation 
from non-spherical SNRs 

SNRs are powerful sources of X-ray radiation, therefore 
X-ray observations of SNR give unique information about 
physical conditions inside the remnants. We calculate here 
the X-ray luminosities in the energy ranges e = hv > 
0.1 keV (L> 01 ) and e > 4.5 keV (L> 4 5 ) as well as the 
spectral index a at e = 5 keV. 



4-.1. Plasma X-ray emissivity under ionization equilibrium 
condition 



We assume the cosmic abundance (Allen |1973| ). The X- 
ray continuum energy emissivity per unit energy interval 
is (in erg cm s" 1 keV -1 ) 



Po(T,e) 



= 1.652 • 10- 23 ^G c (T,e)T e 



-1/2 



exp 



-11.59e 



(14) 



where e > 0.1 is the photon energy in keV, Tq is the 
plasma temperature in 10 6 K, n c is the electron number 
density. The approximation for total Gaunt factor G c as 
the sum of Gaunt factors for free- free Gg, free-bound Gfb 
and two-photon G^ processes was taken from Mewe et 



al. ( P-986D : 

G C (T, e) = 27.83(T 6 + 0.65)" 133 + 0.35e~ a34 T 6 a422 . (15) 

This approximation represents the continuum losses to an 
accuracy of 10 - 20% for T > 3 • 10 6 K and of 30 - 50% 
for T = (0.2 - 3) ■ 10 6 K. 

For calculation of continuum and line emission of 
plasma in the different energy ranges, we have approxi- 



mated the Raymond & Smith (1977) data for plasma emis 

sivity A (in erg cm~ 3 s _1 ) 

for e = 0.1 — 2.4 keV as follows: 

^0.1-2.4/-y>\ _ -lq-21.80ji-0.63 



exp 



-1.40 



(16) 



(accurate to 5 — 30% for logT = 5.3 
for e > 4.5 keV as follows: 



A >4{5 (T) = 10" 



' exp 



-39.50 

rpQ.77 

J 6 



(accurate to 1 — 8% for logT = 6.7 
for e > 0.1 keV as follows: 



A > (T) = A 



0.1-2.4 



(T) + 10 



-24.4rp0.8 



(17) 



(18) 



(accurate to 7 - 35% for logT = 5.5 - 8). 

The total flux F e at photon energy e and luminosity 
L x of the entire SNR can be calculated by integrating over 
the remnant volume V 



F £ = J P c (T,e)dV, 
v 

L x = J A(T)n c nndV. 



The spectral index a is 

d\nF e 
dine 



(19) 



(20) 



(21) 
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Table 1. Luminosity in ranges e > 0.1 keV and e > 4.5 keV (in parenthesis) of SNR expanding in an exponential medium 
Eq. (|ll]) with H = 10 pc (upper lines) in comparison with Sedov SNR (lower lines). The SN explosion energy is E$i = 1. 



log(i yrs ) 










log(rtH) 










-1 




0.5 





0.5 


1 


2.5 


0.056 




34.97 
34.96 


(34.66) 
(34.64) 


35.53 (35.28) 
35.52 (35.27) 


36.12 (35.88) 
36.11 (35.87) 


36.74 (36.42) 
36.73 (36.42) 


3.0 


0.176 


34.65 (34.39) 
34.61 (34.37) 


35.26 
35.23 


(34.93) 
(34.92) 


35.92 (35.42) 
35.89 (35.41) 


36.62 (35.83) 
36.60 (35.83) 


37.35 (36.12) 
37.34 (36.13) 


3.5 


0.555 


35.25 (34.31) 
35.10 (34.33) 


35.95 
35.84 


(34.61) 
(34.63) 


36.68 (34.79) 
36.60 (34.79) 


37.43 (34.83) 
37.38 (34.82) 


38.19 (34.80) 
38.16 (34.78) 


3.75 


0.987 


35.72 (33.94) 
35.47 (33.98) 


36.41 
36.24 


(34.09) 
(34.07) 


37.13 (34.10) 
37.02 (34.05) 


37.86 (34.06) 
37.79 (34.03) 


38.57 (34.11) 
38.54 (34.10) 


4.0 


1.755 


36.23 (33.40) 
35.88 (33.32) 


36.89 
36.66 


(33.41) 
(33.28) 


37.56 (33.37) 
37.42 (33.31) 


38.23 (33.43) 
38.15 (33.40) 




4.25 


3.121 


36.74 (32.86) 
36.29 (32.53) 


37.30 
37.04 


(32.78) 
(32.60) 


37.88 (32.79) 
37.74 (32.72) 






4.5 


5.555 


37.04 (32.54) 
36.65 (31.90) 


37.53 
37.33 


(32.40) 
(32.04) 









4-2. Evolution of the total X-ray emission from aspherical 
SNRs 



It is well-known that the Sedov (1959) solution for SNR 
characteristics in a uniform medium is determined by 
three parameters, e.g., the energy of the explosion E Q = 
10 51 i?5i erg, the initial number density and time t. The 
shape of the integrated spectrum (in particular the spec- 
tral index a) depends only on one parameter (the shock 
temperature T s ) in case of collision ionization equilibrium 
(CIE) and on two (T s and 77 = E a (n° ) 2 ) in case of NEI 
(for fixed abundance) (Hamilton et al. 1983). In order to 
estimate the total luminosity of a Sedov SNR in the NEI 
case, we need a third parameter, i.e. the explosion energy 
E (Hamilton et al. |1983j ). In the Sedov CIE case the lumi- 
nosity depends only on two parameters T s and £ = E^\n\\. 

If the SN explodes in a non-uniform medium, we have 
an additional fourth parameter which characterises the 
non-uniform density distribution. In our cases, it is the 
scale height H for exponential density distribution or r D 
for power-law density distribution. 

It is naturally to expect, that in cases of NSNR con- 
sidered, the total X-ray luminosity L x as well as spectral 
index a will strongly depend on non-uniformity charac- 
teristics, such as scale height H etc. Therefore, we cal- 
culate here extensive grid of models for evolution of X- 
ray radiation from NSNRs which evolve in non-uniform 
ISM and compare it with the Sedov case. Our grid cov- 
ers the range of paramers which correspond to the adia- 
batic stage of SNR evolution, i.e. for SNR radii between 
Ri = (3M cj /47rp°(0)) 1/3 ~ 2.5(n^)~ 1 / 3 pc when the swept 
up mass of ISM becomes equal to the mass of ejecta 
M e j ~ IMq, and Rf when the effective temperature T e f 



Table 2. Spectral index a(5 keV) of NSNR in exponential 
medium Eq. ( [ill ) with H = 10 pc (upper lines) in compari- 
son with Sedov SNR (lower lines). The SN explosion energy is 
Esi = 1. 



log(t yrs ) 








log'K) 






-1 


-0.5 





0.5 


1 


2.5 


0.056 




0.49 


0.58 


0.71 


0.92 








0.49 


0.57 


0.70 


0.91 


3.0 


0.176 


0.75 


0.95 


1.26 


1.74 


2.47 






0.70 


0.91 


1.23 


1.72 


2.45 


3.5 


0.555 


1.95 


2.57 


3.46 


4.52 


5.05 






1.72 


2.45 


3.49 


4.67 


5.06 


3.75 


0.987 


3.06 


3.81 


4.70 


4.89 


4.24 






2.93 


4.10 


5.04 


4.78 


4.22 


4.0 


1.755 


3.80 


4.56 


4.73 


4.07 








4.67 


5.06 


4.45 


4.05 




4.25 


3.121 


4.02 


4.55 


4.01 










4.78 


4.21 


3.94 






4.5 


5.555 


3.60 


3.82 












4.05 


3.85 









corresponds to the maximum value of plasma emissivity 
function T of ~ 1.3T S = 5.2 • 10 5 K (Lozinskaya|1992[). In a 



non-uniform medium, the value of Rf depends on the con- 
crete density distribution in a selected direction (sector). 
So, if the density distribution in a sector is exponential 
n^(r) — 71^(0) • exp(±r/iJ), we may estimate Rf from 
the relation R { ■ exp(±R f /3H) = 25(£ 5 i/n° (O)) 1 / 3 pc. 
For a uniform medium H — 00, we obtain i? max — 
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Table 3. 


The same 


as in Table 1 for 


Eha = 


0.1 and 


range 


e > 0.1 keV. 












log(tyrs) 








Iog(n&) 










-1 


-0.5 





0.5 


1 


2.5 


0.018 


- 


34.12 


34.74 


35.40 


36.10 








34.11 


34.73 


35.39 


36.10 


3.0 


0.056 


33.92 


34.62 


35.35 


36.11 


36.89 






33.89 


34.60 


35.34 


36.11 


36.88 


3.5 


0.176 


34.68 


35.43 


36.19 


36.94 


37.66 






34.60 


35.38 


36.16 


36.92 


37.65 


4.0 


0.555 


35.56 


36.23 


36.86 


37.43 


37.90 






35.42 


36.15 


36.83 


37.43 


37.90 



Table 5. The ratios of SNR characteristics (the volume V, 
swept-up mass M and characteristic temperatures T c h) in 
non-uniform medium (with exponential (E) Eq. ([H]) and power 
law (PL) Eq. @ density distribution) to uniform ones. 



T 


v/v s 


M/M s 


T c h 


/T s 




E 


PL 


E 


PL 


E 


PL 


0.1 


1.00 


1.02 


1.00 


1.00 


1.00 


1.00 


0.5 


1.03 


1.06 


1.00 


0.99 


1.00 


1.01 


1.0 


1.05 


1.10 


1.01 


0.98 


0.99 


1.02 


3.0 


1.14 


1.24 


1.01 


0.95 


0.99 


1.06 


5.0 


1.23 


1.37 


1.02 


0.92 


0.98 


1.09 


10.0 


1.51 


1.64 


1.04 


0.86 


0.96 


1.16 


20.0 


2.53 


2.13 


1.07 


0.75 


0.94 


1.33 



Table 4. The same as in Table | for E 51 



0.1. 



log(tyrs) 








log(n H ) 






-1 


-0.5 





0.5 


1 


2.5 


0.018 




0.71 


0.92 


1.23 


1.72 








0.70 


0.91 


1.23 


1.72 


3.0 


0.056 


1.26 


1.74 


2.47 


3.49 


4.66 






1.23 


1.72 


2.45 


3.49 


4.67 


3.5 


0.176 


3.46 


4.52 


5.05 


4.48 


4.05 






3.49 


4.67 


5.06 


4.45 


4.05 


4.0 


0.555 


4.73 


4.07 


3.87 


3.76 


3.68 






4.45 


4.05 


3.85 


3.74 


3.66 



25^ 3 (^(0))- 



0.3 



pc, which is close to results of other 
In our calculation, the maxi- 



authors (Lozinskaya 1992) 
mal time was estimated from adiabaticity violation in the 
9 = 7r/2 sector. For parameters of ISM considered here, 
maximal radii of SNRs do not exceed a few scale heights. 
Therefore, obtained results are not affected by shock ac- 
celeration, which is important only at considerably greater 
distances. 

We show in Fig. fj] the evolution of X-ray luminosity for 
range e > 0.1 keV and spectral index a at 5 keV for SNR 
in an exponential medium. It is easy to see that both X-ray 
characteristics of NSNR evolve analogously to the Sedov 
case and similarity increases with decreasing surrounding 
medium density gradient (with increasing H ) and decreas- 
ing explosion energy. The results presented in Fig. ^ of 
calculations of the Sedov case for E$x — 1 (Hamilton et 
al. 1983) reveal some differences between X-ray luminosi- 
ties at low temperatures. They may be caused, in part, 
by atomic data and chemical composition differences (our 
model uses the A llen a bundance, whereas in the paper of 
Hamilton et al. ( |l983|) have been used the Meyer abun- 
dance). Our res ults f or Sedov case coinsi de wi th of Leahy 
& Aschenbach ( |l996| ) and Kassim et al. ( 1993|) results for 
model with Allen abundance within 25%. 



We have calculated a grid of NSNR models (Tables |l|- 
^) which confirms that this analogy in evolution is an in- 
trinsic property of NSNR X-ray radiation. In fact, it may 
be seen from Tables that for a wide range of number 
density at the point of explosion n^(0) =0.1 — 10 cm -3 , 
the X-ray luminosity of adiabatic NSNR evolving in a 
medium with strong enough density gradient (H — 10 pc) 
in different energy ranges are not far from the X-ray lu- 
minosity of SNR in a uniform density medium with the 
same initial model parameters. The differences increase 
with age of the SNR and with decreasing n^(0). But even 
for old NSNR (e.g., t — 3 • 10 4 years) in a low density 
medium (n^O) =0.1 cm -3 ), the maximal difference is 
about 60% (for r = 17.5). The differences in spectral in- 
dex of adiabatic NSNR do not exceed 10%. Table |l| shows 
that luminosity in range e > 4.5 keV is close to the Sedov 
case also. These last two facts reveal the similarity of total 
spectra of Sedov SNR and non-spherical ones. 

We suppose that the behevior of integral X-ray char- 
acteristics of NSNR in other cases of smooth continuous 
density distribution of surrounding medium will be simi- 
lar considered above, since X-ray emission mainly depends 
on the emission measure EM ~ n^V ~ A1 2 V~ 1 , but both 
masses of swepted up gas M and volume of NSNR V re- 
main close to those of Sedov SNR (Table |). Moreover, 
the last fact allow us to introduce a characteristic shock 
temperature for NSNR T c h, comparable to T s in the Sedov 
case: 

32tt 7-1 



75 (7 + l) ; 



;(M*)-\ 



(22) 



where T* = T/T m , M* = M/M m with T m = 
a A (2, 7 )~ 1 ^ ga 1 s ^oPo 1 (0)^ m 3 and M m = Po (0)R 3 m . The 
temperature T e f determined from the observed spectrum 
will be conn ected with T c h as in the Sedov case T e f w 
1.3T c h (Itoh 1977 ). This characteristic temperature may 
be used to estimate the parameters of the whole remnant, 
such as mass, volume, etc. 

The results of calculations show that integral X-ray 
characteristics weakly depend on the surrounding medium 
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Fig. 7. Evolution of luminosity L x in range > 0.1 keV and of spectral index a (5 keV) of SNRs in a uniform medium (line 1), 
in an exponential medium Eq. (|ll|) with H = 8 pc (line 2) and H = 10 pc (line 3). In all cases n^Q) = 1 cm -3 . Results of 



numerical calculations of Hamilton et al. (1983) for a uniform medium and E51 — 1 with Meyer abundance are also shown by 
dots. Hereafter 7 = 5/3. 



Table 6. Evolution of SNR luminosity L x and spectral index a(5 keV) in uniform (S) and exponential (E) Eq. (Jill) media. 



logC = log(n H -£ 5 i) 



log(T ch ) 


model 


H 


E51 




-1 













1 




T 


L>^ 


a 5 


T 


i x >U1 




T 


L>^ 


Q5 


8.0 


S 








34.68 


0.81 




35.68 


0.81 




36.68 


0.81 




E 


80 


1 


0.004 


34.68 


0.82 


0.001 


35.68 


0.82 


8e-5 


36.68 


0.81 




E 


10 


0.1 


0.015 


34.68 


0.82 


0.002 


35.68 


0.82 


3e-4 


36.68 


0.82 




E 


10 


1 


0.696 


34.73 


0.88 


0.102 


35.69 


0.83 


0.015 


36.68 


0.82 


7.2 


S 








35.49 


3.02 




36.49 


3.02 




37.49 


3.02 




E 


80 


1 


0.018 


35.50 


3.03 


0.003 


36.50 


3.03 


4e-4 


37.49 


3.02 




E 


10 


0.1 


0.070 


35.51 


3.03 


0.010 


36.50 


3.03 


0.001 


37.50 


3.03 




E 


10 


1 


3.232 


35.75 


3.13 


0.474 


36.56 


3.05 


0.070 


37.51 


3.03 


6.4 


S 








36.56 


4.18 




37.56 


4.18 




38.56 


4.18 




E 


80 


1 


0.083 


36.57 


4.19 


0.012 


37.56 


4.18 


0.002 


38.56 


4.18 




E 


10 


0.1 


0.323 


36.59 


4.20 


0.047 


37.57 


4.18 


0.007 


38.56 


4.18 




E 


10 


1 


15.00 


36.98 


3.73 


2.202 


37.70 


4.38 


0.323 


38.59 


4.20 



density gradient (e.g., H). Therefore, as in the Sedov case, 
a spectral index of equilibrium emission from NSNR de- 
pends approximately only on one parameter T c h and lumi- 
nosity depends approximately only on two parameters: T c h 
and C = E^n^iO). As confirmation of this, we show that 
maximal deviation from this rule for NSNR in exponential 
density distribution reaches few percent for spectral index 
and a few tens of percent for luminosity (Table ||). 



4-3. Discussion 

From Fig. |?] and Tables [I]j6| follows a remarkable fact of 
proximity of total fluxes and spectrum shapes of SNRs in 
uniform and non-uniform media even in case of consid- 
erable anisotropy of adiabatic NSNR. The reason of this 
phenomenon is mutual compensation of emission deficit 
from low density regions of NSNR and enhanced emission 
from high density regions. This explains the strange cir- 
cumstance that integral X-ray characteristics of the ma- 
jority of SNRs with evident asymetry in shape and/or 
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Fig. 8. Distribution of e > 0.1 keV surface brightness S >01 
alon g t he NSNR symmetry axis in uniform and exponential 
Eq. (Jill) media for time moments t = 1000 yrs (line 1), 5000 yrs 
(line 2), 40000 yrs (line 3). The model parameters are E51 = 1, 
rin(0) = 1 cm" 3 , H = 10 pc. 



Fig. 9. Distribution of surface brightness S >01 of the Sedov 
SNR in uniform medium in the range e > 0.1 keV. The model 
parameters are E&x = 1, 71^(0) = 0.2 cm -3 , t = 1000 yrs. 
The SNR characteristics are log(L> 01 , erg/s) = 34.98, 
M = 9.4 M Q , T a = 9.9- 10 7 K. The lines of constant brightness 
are indicated by values of logarithm of flux log(S'). The center 
of explosion hereafter is at the origin of the coordinates. 
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Fie. 10. The same as in Fig. |9| for exponential medium Eq. 
(Ill) with H — 10 pc. The NSNR parameters at this time are 
log(L>°'\ erg/s) = 35.02, M = 9.4 M , T ch = 9.9 ■ 10 7 K. 
The NSNR inclination angle to the plane of the sky equals 0° 
(upper case), 45° (center case) and 90° (lower case). 
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Fig. 11. The same as in Fig. |Io]for time t = 3300 yrs (log(L> 1 ) = 35.70, M = 40.0 M Q , T ch = 2.3 • 10 7 K) and t = 25500 yrs 
(log(L> 1 ) = 37.28, M = 476 M Q , T s = 2.0 ■ 10 6 K). 
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5. Surface - distributed characteristics of NSNR 
X-ray emission 

5.1. Surface brightness in range £ > 0.1 keV 

We analyse here the surface brightness distribution S(x, z) 
(in erg s _1 cm -2 ster -1 ) 

V2 

S{x,z) = i- / A{T)n e n B dy (23) 
yi 

and surface distribution of spectral index 
dine 



a{x,z) = I P c (T, £ )dy, (24) 



yi 



where the integrals are taken along the line of sight inside 
the remnant. 



5.1.1. NSNRs in a medium with flat exponential density 
distribution 

Figs. ||-[ll] demonstrate an evolution of surface bright- 
ness of NSNR in media with exponential density distri- 
bution (p[). At the beginning of SNR evolution, the sur- 
face brightness map is very like the Sedov one (Fig. ^|), 
but with age the differences become more considerable 
when the majority of emission arises from more dense re- 
gions. The brightness contrast S max /S m i n may increase 
with time up to ~ 10 5 (S max and S m - ln are the values of 
both the biggest and the smallest maxima in distribution 
of surface brightness). 

In order to interpret an observation, it is necessary to 
take into account that projection effects also essentially af- 
fect the visible morphology of NSNR. Three cases of pro- 
jection are shown. We can see that projection decreases 
real anisotropy and contrasts. For example, under condi- 
tion of full disclosure of NSNR at the age t = 3300 years 
(Fig. [ll]) S max /Smin = 425 but for inclination angle 45° 
this ratio is S meix / S m i n = 75. 

It is interesting also to compare Fig. || and a bottom 
case of Fig. [lO] (when a visible shape of NSNR is spherical 
as result of projection). One can see that even spherically 
symmetric Sedov-like observational distribution of surface 
brightness does not guarantee an isotropic distribution of 
parameters inside SNR i.e., uniformity of ISM density and 
isotropy of explosion. 



5.1.2. NSNRs in a medium with power-law density distri- 
bution 

Fig. [l^ shows three projections of NSNR in power-law 
density distribution (|l|) with E 51 = 1, n^(0) = 0.2 cm~ 3 , 
t = 1000 yrs. Corresponding NSNR characteristics are 
log(i> 1 ) = 37.03, M = 268 M Q , T ch = 3.5 • 10 6 K. 
It may be compared with values calculated for NSNRs 



Fig. 12. Distribution of surface brightness S l>01 of NSNR 
in medium with power-law density distribution Eq. ( |l3| ) in 
range e > 0.1 keV. Model parameters are: E51 — 1, 
nfi(0) = 0.2 cm" 3 , t = 17800 yrs, r D = -17.5 pc. The NSNR 
inclination angle to the plane of the sky equals 0° (upper case) , 
45° (center case), and 90° (lower case). 
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Fig. 13. The same as in Fig. for range e > 4.5 keV and time t = 3300 yrs (log(L> 4 ' 5 ) = 34.48) and t = 25500 yrs 
(log(L> 45 ) = 32.58). 



which evolve in media with different density distribution, 
to show that different ISM density distributions give simi- 
lar integral (surface-integrated) X-ray characteristics. For 
the same model parameters in case of SNR in uniform 
density, we obtain log(L>° - 1 ) = 36.75, M = 299 M Q , T s = 
3.1 • 10 6 K and for exponential NSNR: log(L> 01 ) = 37.08, 
M = 308 M Q , r c h = 3.0 • 10 6 K. 

As result of projection, the maximum of surface bright- 
ness does not lie close to the edges of NSNR as in shell- 
like SNRs, but creates a compact region inside the visible 
projection. Therefore, the projection effect in case of the 
NSNR elongated predominantly along the line of sight is 
one possible sources of apparent filled-centre SNRs, espe- 
cially when the search for a pulsar has no result. 



5.2. Surface brightness in range e > 4.5 keV 

Different photon energy ranges reveal different sensitivity 
to the non-uniformity of ISM. Fig. |l3| shows the results of 
calculation of surface brightness of NSNR in an exponen- 
tial density distribution (111]) in range e > 4.5 keV. We may 
see that surface brightness contrast in this range is essen- 
tially smaller then in the wider range e > 0.1 keV, where 
line emission dominates. So, for t = 3300 years (Fig. [llj) 
the maximal surface brightness contrast is S max / S m - m = 
425 for e > 0.1 keV and is only Smax/Smin = 7 for 
e > 4.5 keV. A contrast of surface brightness in range 
> 0.1 keV is mainly caused by the contrast of the sur- 
rounding medium density distribution (S oc p 2 ), but in 
range > 4.5 keV, it weakly d epend s on density contrast 
(S cx p 1 ' 2 ) (Hnatyk & Petruk |1996|) . 



5.3. Surface distribution of spectral index 

Plasma in different regions of NSNR is under different con- 
ditions which influence emission. Therefore, spectra from 
different NSNR regions must be different. The spectral in- 
dex distribution for NSNR in an exponential medium is 
shown in Fig. |lj. The contrast of values of spectral in- 
dex increases with time but even for old adiabatic NSNR 
does not exceed a few times. Meantimes, the projection 
effects decrease both the index contrast and anisotropy of 
its distribution. 

5.4- Discussion 

We will summarize the results presented in this section. 

1 . The distribution of surface brightness of NSNR dif- 
fers essentially from a spherical one. 

2. The contrast of surface brightness increases with the 
age of the SNR and may reach a few ten of thousands. 

3. Projection effects hide real contrasts of surface char- 
acteristics of NSNR's emission. Observational morphology 
of NSNR depends essentially on its orientation to the line 
of sight. 

4. A harder X-ray range (e.g., e > 4.5 keV) is less 
sensitive to the influence of surrounding medium non- 
uniformity than a range which includes soft emission. This 
fact may be used for testing ISM non-uniformity. 

5. The surface distribution of spectral index a may also 
be an effective test for NSNR diagnostics. 

6. Conclusions 

Large-scale non-uniformity of ISM density distribution 
in the SN neighborhood essentially affect the evolution 
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Fig. 14. Surface distribution of spectral index a at 5 keV 
for NSNR in exponential medium Eq. ([ll]) with H = 10 pc, 
£51 = 1, n^(0) = 0.2 cm" 3 for time t = 1000 yrs (spec- 
tral index from total SNR a tot (5 keV) = 0.86), t = 3300 yrs 
(a tot (5 keV) = 2.39) and t = 25500 yrs (a tot (5 keV) = 3.95). 



of SNR. The shape of SNR becomes essentially non- 
spherical and distributions of parameters inside the rem- 
nant, as viewed from the center of explosion, are strongly 
anisotropic. We have proposed a new approximate ana- 
lytical method for full hydrodynamical description of 3D 
point-like explosions in non-uniform media with arbitrary 
density distribution i.e., for cases, when well-known self- 
similar Sedov solutions are inacceptable. On the basis of 
it, we carry out the simulation of evolution of 2D non- 
spherical SNRs with special attention to their X-ray radi- 
ation. Since our aim in this paper was to investigate the 
role of density gradients in NSNR evolution, we restricted 
ourself with the ionization equilibrium case in the calcula- 
tion of X-ray radiation. The role of electron conductivity, 
nonequilibrium and nonequipartition effects will be inves- 
tigated elsewhere. 

At first, we have investigated the shape of NSNR and 
brought out a remarkable fact of sphericity of visible shape 
of NSNR, even if the real deviation from sphericity is large. 
When i? max /-Rmm < 2, the observed shape will differ from 
spherical by less than 5%. Visual shape becomes notice- 
ably non-spherical (more than 5% of visual assymetry) 
only for surface density contrast of order 100. Projection 
effects on the sky plane decrease the real shape anisotropy. 

Therefore, we have calculated the equilibrium X- 
ray emission characteristics of NSNR. Again, the total 
(surface-integrated) parameters of X-ray radiation (lumi- 
nosity, spectral index) are close to those in Sedov case 
with the same initial data. This is why many SNRs with 
apparently 2D anisotropic distribution of X-ray surface 
brightness are well described by the Sedov model. 

Contrary to previous cases of weak dependence of 
shape and integral X-ray emission characteristics on den- 
sity gradient, the distribution of X-ray emission charac- 
teristics along the surface of NSNR are very sensitive to 
an initial density distribution arround SN. Moreover, con- 
trast of X-ray surface brightness S caused by density in- 
homogenity is most prominent, i.e., higher than density or 
temperature contrasts along the NSNR surface. For case 
of NSNR in Fig. [lj (t = 3300 yrs), surface brightness con- 
trast is S^^/S^^ 1 — 425, but contrast of shock density 
equals 10 and of shock temperature equals 4 only. 

Therefore, the surface brightness maps give the most 
promising information concerning the physical conditions 
inside and outside NSNR. It is important to note that 
typical contrasts of surface brightness caused by density 
inhomogenity (up to ~ 10 5 ) are considerably higher than 
those caused by nonequilibrium effects. Therefore, the role 
of density gradient is dominant in interpretation of NSNR 
observations. 
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the manuscript. 
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A. Appendix 

A.l. Lagrangian form of hydrodynamical equations and 
shock conditions 

We use the set of hydrodynamical equations for the case of 
one- dimensional adiabatic motion of nonviscous perfect 



gas in Lagrangian form (Klimishin 1984) 



Pt 



P 



p 2 /r\ N 

~\~) Ua 

p° \ aJ 



Nup 



= 0, 



= 0, 



c 2 pt 



0. 



r t — u = 0. 



(Al) 



(A2) 



(A3) 



(A4) 



Here gas pressure P(a,t), its density p(a,t), velocity 
u(a, t) and Eulerian co-ordinate r(a, t) are functions of 
Lagrangian co-ordinate, i.e., initial gas particle position a 
and time t, c — y jP/ p is the adiabatic sound velocity, 7 
is the adiabatic index. Subscripts indicate partial deriva- 
tives with respect to corresponding variables, p° = p°(a) 
is the initial density distribution. 

The continuity equation (A-4) may be written in the 
following form 



-.N 



■ da = p ■ r N ■ dr 



or 



p \r / 



N 



(A5) 



(A6) 



At front of a strong shock with trajectory R = R(t), 
the following conditions are satisfied 



jR 



p* = (p°) s /(l-u;) 



P s = uj{p°) s R 2 



r s = R 



(A7) 
(A8) 
(A9) 
(AlO) 



where R — dR/dt is the shock velocity, superscript "s" 
corresponds to values of parameters at the shock front 
a = R, w = 2/(7 + 1). 



A. 2. Derivatives of functions of parameter distribution at 
the shock front 

Equations (-1) - (-5) and shock conditions (-6) - (-9) allow 
to find the values of arbitrary order partial derivatives of 
hydrodynamical functions at the shock front using the law 
of shock motion R = R{t) (Gaffet |1978[ ). 

To find first derivatives of functions p, P and u at 
front, we use equations (-1) - (-3), written for a = R, 
and add three equations, resulting from differentiation 
of boundary conditions (-6) - (-8) along the shock tra- 
jectory by operator D/Dt — (d/dt) + R(d/da). Solving 
the obtained set of six equations for six unknown partial 
derivatives at the shock front, we obtain expressions for 
<, uf, P a s , P t s , p s a , p\ (Hnatyk & Petruk |996l). 



From equation (-5), we have now 



(All) 



Rr s 



lo{1 - lo) [ZB + N{2 -uj)- ml 



where B = RR/R 2 , m = m(R) — —(dlnp/dlna) s . 

To find second derivatives of hydrodynamical functions 
at the shock front we differentiate equations (-1) - (-3) sep- 
arately with respect to a and i, at the shock front (a=R), 
what gives us six equations for nine unknown derivatives. 
Additional three equations are given by differentiation of 
shock conditions (-6) - (-8) by the operator 



D 2 
Dt 2 



dt 2 



2R 



d 2 



dadt 



R< 



da 2 



da 



(A12) 



Solving the system of nine equations with nine un- 
knowns, we obt ain expressions for u s at , u s aa , p s aa (Hnatyk 
& Petruk 1996D - From (-5) we obtain 



R 2 r s aaa = w(l - lj) [3(7 - 5u;)B 2 + 

+ [(-5a; 2 +4lo + 8)N + {4lo - ll)m] B+ 
+lu(2lj 2 -7lj + 6)N 2 + (lu 2 A)Nm- 
-u(2 - lo)N -{lo- 2)m 2 + (2lo - l)m+ 

+ (2lo - l)m' + (6lo - 4)Q 



(A13) 



where Q = R 2 R^/R 3 , m' = -dm/dlnR = -R ■ 
dm/ dR. 

The time derivative (r s a ) t — 0, and for other ones we 
have 



R 



5 • (Rr s aa ) t = w(l - w) [3(B + Q- 2B 2 ) + m'] (A14) 



R 



5 ■ (R 2 r B aaa ) t = w(l - cj) 

+ [(-5lu 2 +4lu + 8)N + (4lu - ll)m\ ■ 

■(B + Q-2B 2 )- 
-[{4lu - 11) B + (lo 2 +lu - 4)7V + 6(1 - u))m]m' 
-(2w-l)[(2m + l)m' + m"J + 
+2(3u; - 2) [2Q + L- 3BQ] 



3(7 - 5lu) ■ 2B(B + Q- 2B 2 )+ 

(A15) 
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where 

L = = ^Q-2Q + 3BQ 



(A16) 



R 4 R 

Derivatives from B and Q with respect to time are 

B= ^[B + Q-2B 2 ], Q = -[2Q + L-3BQ\. (A17) 
R R 

In the self-similar case [B + Q — 2B 2 ] — 0, 
[2Q + L- 3BQ] = 0. 

According to the approximate formula for shock veloc- 
ity (1) in the common (non-similar) case 



B = k(m)[m - {N + 1)], 



(A18) 



= 2fc 2 (m)[m - (N + l)] 2 - fc(m)[m' + m - (N + 1)], (A19) 



L = (45 - 1)(B + Q - 2B 2 ) + (3B - 2)Q + fc(m)r 



(A20) 



In the work of Hnatyk & Petruk 1996, the rest deriva- 



tives u\ t , p s at , p s tt , P* a , P* t , P t s t are presented. 



r w the r(a) dependence is r/R = C ■ (a/R) x and the 
factor C is connected with central pressure in the following 
way: 



C = C A = 



{3 + N'f 



(7 + l)( 



7 + l\f 



/AT \\ - V(7(A r +l)) 



(A23) 



■R" +L P(0,t)- 



£ 



thus for 7 = 5/3 in case N — (plane shock) Ca = 1.1670, 
and in case N — 2 (spherical shock) Ca = 1.0833. 

In the more general case of anisotropic explosion with 
direction-dependent energy release, E = E(Cl) and shock 
trajectory R — R(Q,t), from (]?]) one obtains generalized 
condition for factor C: 



c(n,t) = 

where E = / E(Q.) dfl. 



/( 7 (JV+i)) (A24) 



j4.5. Approximation of the connection between Lagrangian 
and Eulerian co-ordinates 

The algorithm considered above allows to calculate partial 
derivatives of arbitrary order. We consider the case when 
an expansion of density, pressure and velocity into series 
are restricted by second order at shock front and first in 
the central region. Namely, if at time t the shock position 
is R(t), we approximate a connection between Eulerian co- 
ordinate r(a,t) and the Lagrangian one a in the following 
way 



r(a,t) 
R(t) 



(!)"■< 



l + a-£ + /?.£ 2 +7.£ 3 + <5-£ 4 ),(A21) 



where £ = (R — a) j R, x — (7 — l)/7, and for each sec- 
tor parameters a,/?, 7,(5 are choosen from the condition 
that partial derivatives at shock front r*, r^ a , r^ aa and 
in center of explosion r° correspond to their exact values: 



a — —r„ + x, 



P=^-(Rr s aa -2x-r s a + x(x + l)), 



(A22) 



-3x(l + x)-r s a + x(x + l)(x + 2)), 

6 = C-(l + a + P + j). 

In case of a uniform medium from the Sedov self- 
similar solution, it follows that in the central region at 
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